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Abstract 

As has long been known, it is energetically favorable for massive fermions to 
deform the homogeneous vacuum around them, giving rise to extended bag-like 
objects. We study this phenomenon non-perturbatively in a model field theory, 
the 1 + 1 dimensional Massive Gross-Neveu model, in the large ./V limit. We prove 
that the bags in this model are necessarily time dependent. We calculate their 
masses variationally and demonstrate their stability. We find a non-analytic 
behavior in these masses as we approach the standard massless Gross-Neveu 
model and argue that this behavior is caused by the kink-antikink threshold. 
This work extends our previous work to a non-integrable field theory. 
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A central concept in particle physics states that fundamental particles acquire their 
masses through interactions with vacuum condensates. Thus, a massive particle may 
carve out around itself a spherical region |J or a shell |2| in which the condensate is 
suppressed, thus reducing the effective mass of the particle at the expense of volume 
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and gradient energy associated with the condensate. This picture has interesting 
phenomenological consequences [[l], [|. 

Here we study these effects within the 1 + 1 dimensional massive generalization 
of the Gross-Neveu model |4| (which we will refer to as MGN), 
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describing N self interacting massive Dirac fermions ip a carrying a flavor index a = 
1, . . . , N, which we promptly suppress. As usual, the theory can be rewritten with 
the help of a scalar flavor singlet auxiliary field a(x). Also as usual, we take the large 
N limit holding A = Ng 2 fixed. Integrating out the fermions, we obtain the bare 
effective action 

S[ a ] = -—j d 2 x (a 2 - 2Ma) - iN Tr log (i$ - a) . (2) 

Noting that ^(ip—cr) = —(ip+cr)^, we can rewrite the Tr log(i$— a) as |Tr \og{i$— 
<j)(i0+cr). If a is time independent, this may be further simplified to ~ / f^[(Tr log(/i + — 
uj 2 ) + Tr log(/i_ — uj 2 )] where h± = —d 2 + a 2 ± a' . Clearly, the two Schrodinger op- 
erators h± are isospectral (see Sec. II of 0) and thus we obtain 

S[a] = ~J d 2 x (a 2 - 2Mo) - iN J ^Tr log(/i_ - u 2 ) (3) 

— oo 

In contrast to the standard massless Gross-Neveu model (the GN model), the 
MGN model studied here is not invariant under the Z2 symmetry tp — > 75^, — > —a, 
and the physics is correspondingly quite different. The GN model contains a soliton 
(the so called CCGZ kink [7], ||, in which the a field takes on equal and opposite 
values at x = ±00. The stability of this soliton is obviously guaranteed by topological 
considerations. With any non-zero M the vacuum value of a is unique and the 
CCGZ kink becomes infinitely massive and disappears. If any soliton exists at all, its 
stability has to depend on the energetics of trapping fermions. Also, the GN model is 
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completely integrable, while the MGN model is widely believed to be non-integrable. 
In recent work || |] we have studied integrable models, and one of the purposes of 
this note is to show that it is possible to obtain non-perturbative results even for non- 
integrable models, albeit in the large N limit. 

The vacuum state Setting a to a constant we obtain from (^) the renormalized 
effective potential (per flavor) 
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where \i is a sliding renormalization scale with X(fi) = Ng 2 (/j,) and M([a) the running 
couplings. By equating the coefficient of a 2 in two versions of V, one defined with \i\ 
and the other with /X2, we find immediately that 
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X(fii) Af> 2 

and thus the coupling A is asymptotically free, just as in the GN model. Furthermore, 
by equating the coefficient of a in V we see that the ratio is a renormalization 
group invariant. Thus, M and A have the same scale dependence. 

Without loss of generality we assume that M(/i) > and thus the absolute min- 
imum of (£|), namely, the condensate m = (a), is the positive solution of the gap 
equation 
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Referring to (0), we see that m is the mass of the fermion. Using (|^), we can re- 



write the gap equation as 



which shows manifestly that m, an observable 



A(m) 

physical quantity, is a renormalization group invariant. This equation also implies 
that M(m) = m, which makes sense physically. 

Static space dependent u{x) backgrounds Ideally, we would like to solve the field 



equation 



ss 



Sa(x,t) 



0, a difficult task beyond the capability of field theorists at present. 



A more realistic goal is to restrict ourselves to time-independent a field and to try to 



solve 



ss 

Sa(x) 



0, but even that is difficult since we don't know how to evaluate S for 
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an arbitray time independent but space dependent o~(x). Furthermore, we can show 
(see below) that such a solution does not exist. 

The relevant physics is not difficult to understand. A generic <j(x) will distort the 
fermion vacuum, causing the fermions to back-react on o~(x) so as to minimize their 
energy, and in general a will become time dependent. In our previous work H] we 
have found a necessary condition (albeit generically insufficient) to avoid such a back- 
reaction in 1 + 1 dimensional theories. The condition is physical and easy to state (and 
we will state it in the present context.) Consider the expectation value of the fermionic 
vector current j^(x) in the background specified by a field configuration a. After some 
standard manipulations we could show H that the spatial component of the fermion 
number current will not vanish at spatial infinity, unless the Schrodinger operator h_ 
is reflectionless. Moreover, the fermion number current will run in opposite directions 
at x = ±00. This apparent current non-conservation indicates that a state giving rise 
to a static reflectionful h_ is highly unstable and will immediately try to decay to a 
stable state by emitting fermions. 

We will thus restrict ourselves to only those a configurations which correspond to 
reflectionless Since there are only denumerably infinite number of reflectionless 
Schrodinger operators known, this condition vastly restrict the space of possible a(x). 
Our calculation amounts to a variational calculation in quantum field theory. For any 
given fermion number Nf the energy of the bag or lump we calculate below is an upper 
bound to the true energy. 

A variational calculation of the bag mass This upper bound on the true energy 
cannot be saturated by static a(x) configurations, because as we already mentioned, 
the MGN model does not have static saddle point o~(x) configurations. However, 
it is clear from the discussion above that reflectionless cr(x) configurations are the 
best trial configuration among all static configurations. As usual, the art behind a 
variational calculation consists of a judicious choice of a trial function. 

The energy functional (per flavor) £[cr(:r)] for static a(x) configurations is by 
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definition £ 
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where T is some temporal infrared cutoff. We write 

, I oo 

S[a(x)) = ^Jdx [V(x) - 2Ma{x)\ - J ^ Tr log [-d 2 x + V(x) - 



as 
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(7) 



where V(x) = a 2 (x) — a' (x) . (Here we used / dx a'(x) = by invoking the boundary 

— oo 

conditions a(x) — m.) 

v ' x— >±oo ' 

Out of the denumerably infinite number of reflectionless Schrodinger operators we 
now take the simplest possibility: that h- = —d 2 + a 2 — a' has a single normalizable 
bound state at some positive energy uj\ < m 2 (and thus bound states at ±u)b in 
the Dirac operator.) It is well-known from the annals of quantum mechanics that 
these properties uniquely determine the single parameter family of potentials 

V(x) =m 2 - 2k 2 sech 2 [k(x - x )} (8) 

(up to an overall translation parameter x which we immediately set to zero.) The 
normalized bound state wave function is ipb{x) = sech kx. The bound state energy 
u) 2 is given by k — \Jm 2 — u) 2 , thus suggesting that we trade k immediately for an 
angle | > 6 > such that k = m sin# (and thus cj& = m cos9.) The corresponding 
cr(x) is: 



cr(x) = m + k 



, . 1 m + k\ ( 1, m + k 

tanh [ kx log — tanh [ kx + — log 

4 m — kJ \ 4 m — k 
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With (^|) as a trial configuration, the energy (|7|) becomes an ordinary function 
8 (9). We thus vary with respect to the variational parameter 9 (or equivalently k.) 
The extremum condition on the energy is 
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Here R(x, uj) = (x\(— d 2 + a 2 — a' — uj 2 ) l \x) denotes the resolvent of h-, and can be 
calculated to be HI 



Substituting (0) in (0) we find 
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where I(u b ,m) = J c 



C 2ni (ul~u) 2 ) y/m 2 -w 2 



(the contour C is specified below.) 



Note that the u integral in the first term in dl2|) is logarithmically divergent. This 
UV divergence is taken care of as follows. Let us consider ([12]) at the cutoff scale A. 
Setting a to its vacuum value in SS/Sa(x) = (see ([TJ) below) we have the (bare) 
gap equation 
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Using (|13]) in (O) we see that all reference to A disappears and the extremum condi- 
tion becomes 
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To evaluate the integral I(ub,m), we have to choose the proper contour C, and thus 
we have to invoke our understanding of the physics of fermions. We fill the Dirac 
sea, including the discrete state at — u b , and then put Nf fermions into the state at 
u b . Mathematically, we thus have to let C enclose the cut on the negative u axis and 
then go around the pole at — u b N times and around the pole at u b Nf times. In this 
way, we obtain || § / = (— — u)/m 2 sin2# where we have introduced the "filling 
fraction" u — 

Recalling first order perturbation theory we immediately recognize the matrix- 
element in (|i~4l) as simply du\jde. Putting it all together we find the extremum 
condition 

rlf. \(0 u\ 

sine =0. (15) 
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de =2m 
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(where we have defined the renormalization group invariant ratio 7 = thus fixing 
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6 as a function of the filling fraction 

6 v 

- + 7 tantf = - . (16) 

7T I 

Integrating ( jTjJ ) and using (|T^) we find that the mass M (namely, NS evaluated at 
the extremal point) of our bag or lump is 

M(v,i) 2 . n , 1 + sinfl 

— ^ = - sin0 + 7 log — . (17) 

JMm 7r 1 — smu 

By calculating = — 7rsin#/ (1 + 7T7 sec 2 #) we see that £{v) is a convex function 
and thus satisfies £(^1 + ^2) < £(^1) +£(^2)- Therefore, a lump binding Nv fermions 
is variationally stable against decaying into two lumps with Nu\ and Nv 2 fermions 
respectively (with v = v\ + u 2 < 1.) Thus, the lump binding Nv fermions is the most 
variationally stable static configuration at the sector of fermion number Nv. Note 
that this is true for small as well as for large values of 7. Furthermore, it is clear 
from ([17D that the binding energy (in units of m) per fermion B(u,j) = 1 — 
increases with u, and does not saturate as in nuclear physics. This is characteristic 
of the bag picture, in which each additional particle digs a deeper hole in the vacuum 
condensate ("the mattress effect"). To demonstrate these facts we present in Fig. (1) 
a numerical computation of the binding energy per fermion at a particular 7. 

binding per fermion 




Figure 1: The binding energy per fermion £?(i/, / y) at 7 = 0.1 
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Let us consider the most stable bag, namely the bag at v — 1. In the small 7 
limit, 

^)~!- 7 log^ + 0( 7 i). (18) 
Nm 7r 4 

Thus M(l, 7) is non-analytic at 7 = 0, i.e., at the GN point. Note that M(l, 0) = ^ 
is the kink-antikink threshold of the GN model [[J. It appears that the logarithmic 
singularity in M(l,7) is associated with the enhanced Z 2 symmetry at 7 = 0. To 
further argue in this direction, we find that as soon as v decreases from 1, which puts 
us below the kink- antikink threshold, M(u, 7) = ^cos 71 ^ ^ + [log ( ^ ) — 

- ^)sin^% + 0(7 2 ) = I - 2 7 log ^ + 0((1 - ^) 2 , 7 2 , (1 -^) 2 7) • 
Thus, the logarithmic singularity 7log7 as 7 — >• is replaced by a logarithmic singu- 
larity 7log(l — v) as — ► 1. This means that the kink-antikink is indeed the source of 
this singular behavior. It would be interesting to address this issue in the framework 
of an appropriate effective action for bags. 

We now turn to the large 7 limit, which may be attained by making the four-fermi 
interactions weak. The theory should then describe quasi-free heavy fermions of mass 
m. We thus expect that the binding energy of bags will tend to zero as 7 — ► 00. This 
is indeed the case, and we find M ^' 7 - > ~ v — 7^ (hi) + C(7~ 3 ). The u 3 behavior is once 
again a manifestation of the mattress effect. We present the results of the numerical 
computation of the binding energy per fermion at maximal filling B(l, 7) in Fig. (2). 
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Static bags do not exist The extremum condition SS/Sa(x) = reads 

a(x) — M r . . „ . °r du> . 1 . 

i ~z = 2a x + d x / — x 2 5 x . 19 

A J Z7r —a; + cH — a' — c^r 

-oo x 

A static saddle point configuration o~(x) is necessarily reflectionless. Concentrating 
on the case of a single bound state and substituting (|TT|) into (|19D and using (|I3D , we 
find 

(cr-m) = (2<r + a x ) (m 2 -a 2 + ( T , )(-^(^ ) m )/ 4 )- (20) 



m\(n) 

Further substituting o~(x) as given in (|9|) into this equation, we verify easily that there 
is no combination n(m, M, A) for which fl20|) is satisfied. Thus ([]]) does not have a 
static reflectionless a(x) saddle point with a single bound state. (The possibility that 
it has a static reflectionless saddle point with more than a single bound state seems 
highly unlikely, but we have not ruled it out rigorously.) We thus conclude that 
the MGN model does not have any static o~(x) configurations. This is in contrast 
to the GN model, which has, as we have already mentioned, static topological o~(x) 
configurations, all of them are reflectionless, with a single bound state f| [7|, ||, or 
more 
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The time dependent bags in the MGN model may be seen as continuous vibrating 
deformations of the static non-topological bags of the GN model (at least for small 
values of 7.) Indeed, the bag configuration in the GN model with one bound state 
is also described by the a field given in (^). Dashen et al [H] showed long ago that 
Kj and m are quantized according to ( |T6| ) and (|T7D, with 7 set to zero || [7]. These 
quantization conditions, as well as the gap equation (|5p, are continuous at 7 = 0. 
Thus, this configuration should not change abruptly as we turn M on, but it cannot 
remain static. This means that as we switch on 7, o~(x) will start vibrating around 
the static profile (^). The vibration amplitude and frequency of these objects must be 
continuous functions of 7 that vanish as 7 — > 0. It would be interesting to determine 
whether and how the lack of static o~(x) saddle points in the massive GN model is 
related to the common lore that turning M on destroys the complete integrability of 
the GN model. 
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The fermion current operator and the bosonized theory By using the methods of 
Sec. 2 of we easily obtain the expectation value of the conserved fermion current 
= tpj^ifj in the background of a static extremal o~(x) configuration trapping Nv 
fermions. The spatial component is identically zero and the fermion density (a j° a) 
is 

p(x) = ^ (a 2 - m 2 ) , (21) 

oo 

which has the correct normalization J dx p(x) = Nv, as can be seen from (^). This 

— oo 

means that, in the bosonized theory, the flavor singlet boson <p develops a spatially 
varying profile which follows the profiles of cr(x) according to 

W v , 2 2 



dMx) = ^~ - (a 2 -m 2 ). (22) 
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